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Abstract 

We present a method to analyze the suitability of particular photonic 
Q-l cavity designs for information exchange between arbitrary superposition 

states of a quantum emitter and the near-field photonic cavity mode. 
As an illustrative example, we consider whether quantum dot emitters 
embedded in "L3" and "HI" photonic crystal cavities are able to transfer 
r/3 ■ a spin superposition state to a confined photonic superposition state for 

use in quantum information transfer. Using an established dyadic Green's 
. function (DGF) analysis, we describe methods to calculate coupling to 

arbitrary quantum emitter positions and orientations using the modified 
local density of states (LDOS) calculated using numerical finite-difference 
time-domain (FDTD) simulations. We find that while superposition states 
are not supported in L3 cavities, the double degeneracy of the HI cavities 
supports superposition states of the two orthogonal modes that may be 
described as states on a Poincare-like sphere. Methods are developed to 
comprehensively analyze the confined superposition state generated from 
an arbitrary emitter position and emitter dipole orientation. 

1 Introduction 

Quantum emitters such as self- assembled III-V quantum dots (QDs) have long 
been proposed as a key component in quantum computing architectures. Many 
of the specific applications of these QDs rely on transferring information between 
the polarization of a photon and the carrier spins in exciton states, that are made 
possible due to strict atomic-like selection rules. The selection rules governing 
interaction of photons with the QD ground state transition are well-defined. 
Due to angular momentum conservation, a circularly polarized photon with 
right (left) helicity will only interact with an electron in the | f) spin ±1/2 
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state and a hole in the | ft), | 4J-) ±3/2 spin state, respectively, giving angular 
momentum change Al — ±1. 

This selection rule is the basis of schemes to use QDs to store and manipu- 
late photon polarization states using carrier spins. Of several widely-proposed 
components we discuss one in this paper. A QD spin-photon interface [T] relies 
on the optical transitions between a resident electron spin and a charged exciton 
(or "trion" ) state. Indeed, any superposition of states on the electron spin Bloch 
sphere, a\ f) + P\ I) has a one-to-one correspondence with a photon polariza- 
tion state on the Poincare sphere, a\L) + (3\R), where a and (3 are complex, and 
contain the crucial phase information for quantum information. The selection 
rules have been shown to be accurate for light interacting with the QD in the 
paraxial approximation, and in a vacuum with uniform local density of states 
(LDOS). 

We note also that the same approach may be used for another application: 
a QD entangled photon pair source produced via the cascaded decay of the 
QD biexciton-exciton [2 . This source produces two photons entangled in linear 
polarization, whose fidelity is dependent on the x and y exciton decay paths 
being indistinguishable, adding an extra complication to the analysis [3]. 

In order to increase the light-matter interaction strength for the high demands 
of quantum information technology the QDs must be incorporated into a photonic 
structure. Along with micropillar cavities [2j, photonic crystal cavities (PhCCs) 
have become a popular means to increase light-matter interaction. Their high 
quality (Q)-factors and low modal volumes have already allowed the strong 
coupling regime to be reached [5]. 

The use of these PhCCs for any applications which rely on the spin selec- 
tion rules is problematic, however. The light matter interaction is increased 
by strongly modifying the LDOS in the cavity. The cavity produces a highly 
anisotropic LDOS. This means that a QD interacting with the LDOS couple 
to one polarization state more strongly than the other. Moreover, the phase 
information contained in the spin state a\ t) + @\ i) may be lost if the coupling 
to the | f) and | |) differs. 

The QD interacts with the intrinsic modes of the cavity, which also have a 
polarization-dependent spatial variation. The position of the QD is also crucial: 
the E-ficld spatial distribution of each polarization component differs, even for 
the case where the cavity produces two modes degenerate in energy, as we 
demonstrate. Finally, the trapped photon then eventually either leaks into the 
far-field outside the structure, or is transmitted down a waveguide. The polar- 
ization of the light collected and measured therefore depends on the difference 
in photon lifetime and angular emission for each polarization component. 

Thus it is clear that coupling photons to spin polarization states in cavities 
is more complex than in the paraxial approximation, where a simple one-to-one 
correspondence between the Bloch and Poincare spheres exists. In this paper 
we discuss a methodology to analyze PhCC designs for their suitability for spin- 
dependent quantum information processes in the weak-coupling regime, taking 
as examples well-known "L3" and "HI" cavity designs. 

We will discuss that in order to use photonic crystal cavities to transfer spin- 
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polarization information with high fidelity, the cavity must meet two criteria. 
The first, Criterion 1, is that the cavity should have doubly-degenerate modes 
(or at least spectrally overlapping modes). The doubly-degenerate modes are 
orthogonal states, allowing an arbitrary superposition state to be generated. 
This allows one to build up a Poincare-like sphere, such that one is able to 
excite the two modes (equivalent to the "H" and "V" linear polarization states 
in the paraxial approximation) with arbitrary intensity and phase. In fact, 
this criterion has been addressed experimentally by applying various methods 
to tune fabrication imperfections in HI cavities to produce degenerate modes 
[51(7] ■ Criterion 2 is that the modes, which have a complex structure in the x-y 
plane, should possess regions where the LDOS is equal for both modes. This 
implies that the QD would couple equally to both, such that there is a one-to- 
one correspondence between the QD Bloch sphere states, and the Poincare-like 
states of the near-field. It is this criterion that is almost completely unstudied 

®- 

We will demonstrate methods to determine whether the criteria above are 
met, and use as examples the well-known L3 and HI PhCC cavities. We will 
show that polarization-degenerate coupling to an L3 cavity mode is not possible, 
whilst a small region in the exact center of the HI cavity shows polarization- 
independent emission. 

The layout of this paper is as follows: In Section [3J we will discuss both 
classical and semiclassical techniques to find an expression for a position and 
orientation dependent response of the QD to the cavity field. In Section [3l we 
will briefly describe the properties, such as dimensions, modes, wavelengths and 
Q factors, of the HI and L3 cavities we will be investigating. We will then apply 
this technique in Section [4] and show the differences between the two types of 
cavities. Note that in this paper we deal with the immediate near-field inside 
the cavity only. 

2 Modelling the light-matter interaction 

Modification of the spontaneous emission lifetime of an excited state of a quantum 
emitter, often known as the Pureed effect [8], has different interpretations in 
classical and quantum physics, both of which are worth considering here. From 
a classical point of view, one may consider the system as a current source dissi- 
pating energy in the form of electromagnetic waves. The waves are scattered 
by the surrounding structure and arrive back at the source, stimulating it to 
dissipate more energy. From a quantum physical point of view, spontaneous 
emission is caused by stimulated emission from the vacuum field. By placing an 
atom in an optical cavity, where the vacuum field is altered, the spontaneous 
emission of the atom can either be enhanced or inhibited. 

In this section we will first show different well-established methods, both 
classical and semiclassical, to describe the basics of the light-matter interaction. 
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2.1 Theory 

The most common approach to this problem is a semi-classical method using 
Fermi's golden rule [5HTT]. Fermi's golden rule uses time-dependent perturba- 
tion theory to find an expression for the decay rate from an initial state \i) to a 
final state |/): 

1 =-5\(f\H i \i)\ 2 p(u;) (1) 
h 

where if ; = /x- £ (7*0) is the interaction Hamiltonian in the dipole approximation, 
£(ro) the vacuum field at the source location, /x the dipole moment and p(oS) the 
density of available optical states. The well-known Purcell factor is then given 
by the ratio of the decay rate in the cavity to that in free space (70) plus the 
enhancement due to leaky modes of the photonic crystal and cavity (F pc ) [10] : 

Fv= 7 _ 3Q(A c /n)3 1 

where Q is the quality factor of the cavity, A c the resonant cavity wavelength 
in free space, n the refractive index of the medium, V m the modal volume 
as given by Eq. T] the spatial mismatch term given by i^nffi^j where 
E(ro) is the local electric field at the source location and E max is the maximum 
field strength in the cavity. We point out that E(r) is a classical field that 
obeys Maxwell's equations. The Lorentzian term pp4QW^j7^ describes the 
spectral mismatch of the dipole with the cavity resonance. When the dipole is 
perfectly aligned with an anti-node of the cavity field and spectrally with the 
cavity mode, Eq. ([2]) reduces to: 

Pmax - 4n 2 V m 1 ' 

The modal volume, V m , is determined using [T2] : 

1 _ R f 1 \ ((/cl/c)) 

V m ~ \v m r Vm - e(r c )mr c ) W 

where f c is the cavity mode calculated from a wave equation with outgoing wave 
boundary conditions and r c denotes the location of the cavity field antinode. 
The norm ((f x \fy)) is defined as: 



«/x|/„» = Jim f e(r)f x (r) ■ f y {r)dr + 1^— f f x (r) ■ f y (r)dr = 6 x , y 

V ^°°JV Ux+UJyJgy 

(5) 

where uj x . y are complex frequencies due to the boundary conditions used and 
SV is the border of the volume V. It should be noted that Fermi's golden rule is 
only valid in the case that the spectral linewidth of the dipole is much smaller 
than the cavity linewidth pj] and in the weak coupling regime. Although these 
two conditions are met in the scope of this work, we will not use the framework 
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of Fermi's golden rule for the analysis in this paper. The complete analysis can 
however be performed within this framework (see analysis for Fig. [6]). 

Alternatively, we may calculate the Purcell enhancement from the ratio of 
energy dissipation by an infinitely small current source in the photonic structure 
to that in a bulk material [5]. The rate of energy dissipation by a current source, 
dW/dT, assuming a harmonic time dependence for both the current source and 
electric field, is given by: 



dW 
~df 



±J Re{j*-E}dr (6) 



where V s is the source volume, j the source current and E the electric field 
vector. Following [9], we can write the source current as j(r) = iu}/i5[r — r ], 
where \x is the dipole moment of the source and r*o is the center of the charge 
distribution of the dipole. Using this we can rewrite Eq. ([6]) as: 

f = -|W™ } (7) 

Alternatively, we may write the electric field created by the source in terms of 
its dyadic Green's function (DGF) [S]: 

E(r) = uj 2 ^ G(r, r ; w)n (8) 

where fj, is the relative magnetic permeability, fiQ the vacuum permeability and 
G(r,ro;uj) the DGF. It is defined as: 

uj 2 

V x V x G(r , r Q ; uj) - — e(r)G(r, r ; uj) = IS{r - r ) (9) 
&■ 

where c is the vacuum speed of light, e the relative electric permittivity and / is 
the unit dyad. Here, G both includes transverse and longitudinal contributions 
e.g. G — G T + G L . Using Eqs. (|7|) and © we can rewrite the energy dissipation 
of the current source as: 

where ft, is the normalized dipole moment and eo the vacuum permittivity. Note 
that here we slightly differ from [5] in order to take complex dipoles into account. 
We can then use Eq. (fTU)) to calculate the energy dissipation of the source. The 
ratio between the power dissipation of the source in the cavity and in a bulk 
material is then given by: 

F P (ro) = ^^"'^^ (ID 
lm{fx*- G {r (h r ;uj)- fj,} 

where Im {Go(tq, Tq; uj)} = ui^Je/Qirc is the DGF in a bulk material. The Purcell 
factor is the normalized equivalent of the semi-classical LDOS, where the LDOS 
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is also proportional to the imaginary part of the DGF fT5[. We can now intro- 
duce the transverse DGF for an arbitrary cavity structure at the dipole location 
following [H[T5]: 

G- { r Q ,r^) = i^^l (12) 

where w c is the cavity resonant frequency and e c is the normalized cavity mode 
function given by: 

e c (r) = E c {r)/y/v^mayi(y/~eE c (r)) (13) 

Because we only consider the transverse part of the DGF here, we do not take 
coupling of the QD with leaky modes into account. Combining Eqs. (TTT1) . (|12l) 
and the GDF in bulk material finally results in: 

„ , > 6ttc 3 f „ e c (r ) ® e*(r ) J 

Fp{r ) = — W M . . (14) 



2.2 Modelling of the polarization-dependent QD-cavity inter- 
action 

Equations ([B1ll4p give a general description of how dipole emission is modified in 
a cavity system due to the change in light-matter interaction. These equations 
do not directly account for the time-dependent decay dynamics of QDs but 
merely take a dipole moment into account. We are thus making the assumption 
that the QD in the cavity may be modeled as a superposition of dipole emitters 
with particular orientations in the x-y plane. Taking the trion transition in a 
QD as an example, a trion in the superposition state ^/l/2(\ f, 4; 1*1") ± | t> U -H-)) 
emits an H(V) photon, leaving an electron in superposition state -\/l/2(| f) ± | 4- 
)), corresponds to a dipole emitter oriented along the x(y)-axis, respectively. 
Indeed, one may describe any arbitrary transition on the spin Bloch sphere 
as two dipole emitters oscillating with a given relative phase. For example, a 
spin I f) (coupling to \L) circular polarization) may be modelled by using two 
dipole emitters that have a relative phase of +tt/2. However, when the strong 
coupling regime or QD-QD interactions are of interest, the time dependent decay 
dynamics of the QD have to be taken into account [THIHZ] ■ 

To calculate how a dipole emitter with a particular orientation and position 
in the cavity will couple to the cavity modes, one may use the Green's function 
method described in [IS]. As we will demonstrate, using the Green's function 
method is more straightforward than using Eq. ([2]), as the situation is not 
trivial for situations such as degenerate cavities. The Green's function method 
describes, for a particular QD orientation and direction, the complete coupling 
to the electromagnetic environment. The approach in [IE] also takes strong 
coupling, coupling of several emitters and several mode frequencies into account. 
In the following section we show, using a similar Green's function method, 
that for a spin-photon interface, one must split the Green's function into two 
degenerate, orthogonal states that have a one-to-one correspondence with the 
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QD spin-transition. We show that an (x,y) "map" of the orthogonal states is 
beneficial in identifying optimal positions in the cavity, and whether cavities are 
at all suitable as a spin-photon interface. 

3 Properties of the investigated cavities 

There are several types of optical microcavities used for coupling QDs to a 
controlled photonic environment. Several designs have demonstrated strong 
coupling, such as circular [TH] and elliptical [3D] micropillars, photonic crys- 
tals [5JIH], and microdisks [32]. The methods we apply below (identifying the 
photonic superposition states, and the correspondence between the spin and 
superposition states) will apply to any lossless dielectric cavity system. We 
choose to study photonic crystals in detail here. This is because these repre- 
sent the most complex case: a circular micropillar or microdisk has a radially 
symmetric LDOS which is relatively easy to analyze in terms of its polarization. 
Even upon lowering the symmetry (for example in an elliptical pillar) the LDOS 
varies more slowly spatially (nevertheless, the polarization-dependent coupling 
of a QD to these now-polarized modes is highly non-trivial [20]). Photonic 
crystal cavities, on the other hand, do not show rotational symmetry of their 
degenerate modes, and show a strong spatial dependence in the polarization 
of the LDOS, as we demonstrate. The lack of radial symmetry makes iden- 
tification of suitable superposition states non-trivial, as these states also do 
not show rotational symmetry. More recent focus on plasmonic-QD structures 
(e.g. [2"3"H2"5] ) suggests that ultrasmall (<< (n\) 3 ) metal plasmonic cavities and 
antenna structures may also be used, again with strong coupling theoretically 
possible. However, these metal structures are dispersive. This means that the 
techniques used here do not directly translate to plasmonics, or lossy systems, 
as the mode expansion is not as straightforward as described here [26]. However, 
we hope to address this in future work. 

By way of demonstration, we will investigate two PhCC designs, both shown 
in Fig. [1] The L3 cavity in Fig. Ufa) consists of a line of three missing holes. 
The dimensions of the L3 cavity we used are based on those presented in [37J , 
where the end holes of the cavity are displaced to obtain a smoother transition 
between the defect and photonic crystal, resulting in a higher Q-factor. The L3 
cavity has the following parameters: a hole radius of r = 0.28a, where a is the 
lattice constant, a slab thickness of d = 0.58a and an end hole displacement of 
s = 0.15a. This cavity design is well-studied, and was the first PhCC design 
to produce strong QD-cavity coupling 5 . This cavity design, however, is not 
known to show any degenerate mode and, experimentally, all modes measured 
from this cavity are strongly polarizing |27[|2"8] . We therefore already do not 
expect that this design is suitable for polarization-degenerate coupling. 

The HI cavity, based on the design in [29 , on the other hand, shows six 
fold rotational symmetry, Cq v , giving rise to an energy degenerate fundamental 
mode [30 . The exact design is shown in Fig. [ljb). The structure has a hole 
radius of r = 0.35a, a slab thickness of d = 0.5a, an end hole displacement of 
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s = 0.10a, the radius of the end holes is r& — 0.25a. All cavities are created in a 
GaAs slab membrane with a refractive index of n = 3.4. Because the photonic 
crystal is realized in a slab, only a bandgap for TE-modes is present and thus 
only TE-modes are confined by the PhCCs, meaning the modes only have H z , 
E x and E y field components [3Tj . 
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Figure 1: Schematic representation of the (a) L3 and (b) HI photonic crystal cavities. The 
teal circle indicates the altered holes surrounding the cavity. The dotted circles indicate their 
unaltered size and position. 



We simulate the structures with the finite difference time domain (FDTD) 
method [32] . using a modified version of an opensource software package [33] . 
A lattice of 13a x 9a and 11a x 11a for the L3 and HI cavity respectively is 
used, we chose a resolution of 90 grid points per lattice period. The first step 
is to find the frequency of any high-Q modes with which the QD interacts. 
This is performed by producing a spectrally wide current pulse, a normalized 
width of Aw c /td c = 0.1 is used, allowing a wide range of modes to be excited 
simultaneously. The modes and their Q-factors are identified using Fourier 
transform and a harmonic inversion technique [34| . We then proceed to extract 
the amplitude and phase information of the mode by outputting time harmonic 
fields at the center of the mode resonance, using a run time Fourier transform. 

An overview of the L3 cavity modes and their properties can be found in [27) 
and in 29 for the HI cavity. In this paper we focus on the fundamental [ll] 
mode in the L3 cavity (as shown in Fig. ^ which is of most interest due to 
its high Q-factor and the dipole mode in the HI cavity (as shown in Fig. [3]) 
which is of interest due to its double degeneracy. The Q-factors found in our 
simulations, 7600 for the fundamental L3 mode and 1300 for the HI dipole mode, 
are relatively low compared to the 47,000 and 15,000 mentioned in [27j[29] for 
the L3 and HI fundamental modes. This is due to the limited size of the PhC. 
By adding more periods to the photonic crystal an increase in the Q-factor can 
be obtained. However, considering the large resolution of 90 grid points per 
lattice period used to simulate these structures, the computational cell needed 
to be limited to keep computation feasible. After testing with larger structures 
in our simulations, we are confident that the reduced size of the PhC does not 
noticeably affect the polarization properties we study in this paper. 

Figure [2] shows the cavity electric and magnetic field distributions of the 
L3 cavity, generated by placing a dipole source in the center of the cavity. 
We extract the field distributions at the peak wavelength of the fundamental 
mode and show that it agrees well with calculations obtained using guided mode 



8 



expansion techniques |27] . 




Figure 2: Cavity field distributions of the L3 mode (a) \H&\, (b) \E X \ and (c) \E y 



We notice from Fig.[2jb,c) that the fundamental mode contains both x and y 
components. However, their spatial profile differs greatly, which strongly affects 
the coupling of the QD to the mode. For example a QD placed directly in 
the center will only couple if it has an y component to its dipole, as the E x 
component at the centre is zero. Likewise, there are regions of the cavity where 
the QD dipole should have an x component to couple to the cavity, as the E y 
component is zero. It already becomes clear here that polarization information 
is not preserved when a QD couples to the cavity. The orientation of the QD's 
effective dipole will determine the rate of coupling to the cavity, as given by 
Eq. (H3J). However, all photons emitted by the QD will emit into the mode, 
which has its own, independent polarization profile, and none of the original 
spin/polarization information remains. 

Thus it would already appear that the L3 cavity is unpromising for spin- 
based QD applications. More promising is the HI cavity. Unlike the L3 
cavity, the HI cavity has two sets of doubly degenerate modes, the "dipole" 
and "quadrupole" modes, as well as two non degenerate modes, the "hexapole" 
and "monopole" modes [29]. Figure [3ja-d) shows the \E X \ and \E y \ components 
of the two dipole modes, as calculated by placing an x and y oriented dipole in 
the center, for the mode in Fig. (3{a,b) and Fig.[3Ic,d) respectively. 

All these modes, just as those in the L3 cavity, are eigenmodes of the PhCC. 
Eigenmodes have the property that a linear combination of the two modes is 
also a valid mode. As we will see, this makes these modes much more viable for 
use in spin-based QD applications. Geometry effects in environment-dependent 
coupling is discussed in detail in [13], where the spontaneous emission rate of 
a dipole in the presence of an arbitrarily-shaped nano-photonic environment 
is analyzed. Counterintuitively, there is no simple relationship between the 
symmetry of the photonic environment and the directional dependence of the 
Purcell enhancement. For example, as the HI cavity has a six-fold degeneracy, 
one might expect a six-fold degeneracy in the DGF, and therefore the Purcell 
enhancement. However, because the DGF may always be diagonalized, it will 
have three eigenvalues (g a , gt>, <7c)> that correspond to three eigenmodes of the 
system. In our system however, because we are dealing with TE-modes only, 
the third eigenvalue will always be zero and can be neglected. When considering 
Purcell enhancement in the cavity, these eigenmodes correspond to two principle 
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dipole orientations with maximum and minimum enhancement, T max ,T m i n . 
The emission rate does not inherit the symmetry of the cavity system. Note 
also that the DGF is dependent both on the nanophotonic environment and the 
position of the dipole emitter. 

In our system, we wish to find cavity designs which have QD positions where 
^max = Tmiri- From the analysis above, this would mean that the two eigen- 
values (g a ,9b) are equal. From Eq. (fT2t it can be easily shown that the eigen- 
values of the GDF can only be equal in the case that two orthogonal modes 
are present and overlap spectrally, the total GDF is then simply the sum of the 
GDF of the orthogonal modes, e.g: 



where e cn is e c of the n cavity mode. The second requirement, finding optimal 
QD positions within the x, y plane is more difficult, and we will discuss this in 
Section HJ However, the first requirement that the eigenvalues are the same 
naturally implies that two orthogonal and degenerate modes must be present. 
We therefore need to first determine whether this is the case for the particular 
cavity in question. 

The dipole mode, which is the only one we consider in this paper, may be 
described in many different orthogonal basis vector fields. Two vector fields are 
said to be orthogonal if their inner product is zero. The inner product between 
two vector fields is given by: 



With this knowledge we can construct different sets of basis vector fields that 
describe the dipole mode of the HI cavity. We choose the HI dipole mode as it 
has been studied fairly widely previously 3,30. These studies have shown that 
two degenerate dipole modes exist in the HI cavities. These may be identified 
by FDTD simulations using an x and y polarized dipole emitter in the center 
of the cavity. Figure EJa-d) show the \E*(x, y)\, \E*(x,y)\ field distributions 
resulting from an x-dipole, and the \E]j!(x, y)\ and \E^(x,y)\ field distributions 
resulting from a y-dipole. To avoid confusion with the field components, we 
have named the different orientations of the mode with a Greek letter. For 
instance the x-mode, noted in terms of its magnetic vector field as H x , is the 
mode that can be excited with a x-dipole in the center of the cavity, while the 
V'-mode is excited with a y-dipole in the center of the cavity. 

Just as diagonally and anti-diagonally-polarized photon states may be constructed 
as the superposition of an x and y-polarized photon, so the "5" (diagonal) and 
"a" (antidiagonal) modes may be constructed with the superposition of the \- 
mode and ^-mode eg, H s / a (x,y) = (H x (x,y) ± H^(x,y)) /y/2. These modes 
are excited by diagonal and anti-diagonal dipoles in the center of the cavity. 
Note that one does not need to re-run the FDTD simultion, as it is iden- 
tical to adding the vector fields. Similarly, the p and A modes are formed 




(15) 




(16) 
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from a phase delayed superposition of the %- and ?/>-modes: H p / X {x,y) — 
(H*(x, y) ± iff* (re, y)) /V2 and can be excited using a right or left circular 
polarizing dipole. One may form a Poincare-like sphere on which all possible 
states of the cavity mode can be mapped. The modes denoted by the Greek 
letters x, ?/>, S, a, p, A, are near-field equivalents of the far-field paraxial x, y, d, a, R, L 
polarization states on the conventional Poincare sphere, respectively. Figure[3Ja) 
shows a Poincare-like sphere, which shows the \H z (x, y)\ field distribution for 
the bases discussed above. Just as the Poincare sphere's surface represents 
the complete set of possible photon polarizations in the far-field, this sphere 
represents the equivalent for the near-field in the cavity and has a one-to-one 
correspondence to the spin Bloch sphere. The ability to form a complete set of 
states on the Poincare-like sphere is crucial for quantum information applica- 
tions. It must be possible to create a set of photonic cavity states that allow a 
full set of superposition states to be contained with full fidelity in the cavity. 




Figure 3: Cavity field distributions for the HI cavity (a) \E X \ and (b) \E y \ field distributions 
of the x dipole mode and (c) \E X \ and (d) \E y \ field distributions of the tp dipole mode, (e) 
Poincare-like sphere with different orientations of the \H Z \ component of the HI dipole mode. 



4 Enhancement in L3 and HI cavities 

So far we have examined L3 and HI cavity modes and shown that for the 
energy-degenerate HI cavity, it is possible to generate a full set of photonic 
modes that would represent superposition states. However, the position of the 
QD determines how much of each component is excited. We examine this in 
this section. Using the results from Section [5J we are now able to look at the 
effects on the Purcell enhancements of an emitter in a PhCC. We will start with 
the [II] mode of the L3 cavity followed by the dipole mode of the HI cavity. 
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4.1 Enhancement in L3 cavities 



From Eq. (TUT) we know that the enhancement is proportional to the square of 
the electric field. The cavity field of the [II] mode can be seen in Fig. [2Jb-c). 
We now plot the enhancement for both x and y oriented dipoles (Fp x (x,y)) 
and {Fp y {x,y)) as a function of position, as shown in Fig. UJa-b) assuming 
that the QD transition is on resonance with the cavity mode. We can also plot 
the enhancement for diagonal (c) and anti-diagonal (d) right handed (e) or left 
handed (f) oriented dipole sources by rotating the E x (x,y) and E y (x,y) fields 
to their Ed(x,y), E a (x,y), E r (x,y) and Ei(x,y) counterparts, i.e. E d / a (x,y) = 
(E x (x,y)±E y (x,y))/V2 and E r/l (x,y) = (E x (x,y)±iE y (x,y))/V2 and calcu- 
lating the GDF using Eq. ([14]) in these bases or by simply rotating the dipole 
moment. From these six graphs we can already clearly tell that the cavity-dipole 
coupling is highly direction dependent for linear sources. For instance at the 
peak of the enhancement Eq. ([3]) tells us the Purcell factor is about 750. This 
point is exactly in the center of the cavity where there is a node in the E x 
field component and an antinode at the E y component. This means that the x 
component of a dipole will not couple to the cavity mode and, depending on the 
F pc factor in Eq. ([2]), emission might be completely suppressed. On the other 
hand, for a circular polarized source the Purcell enhancements are identical for 
both left and right circular polarization. 

To gain a better understanding of the polarizing properties of the cavity, we 
plot the degree of polarization (DOP): 



DOP hJ {x,y) = ZTA~VJ)J..{ ( 17 ) 



Fpijx^y) - Fp 3 {x,y) 
Fp t {x,y) + Fpj(x,y) 

where i, j are orthogonal bases {x, y; d, a; r, I), in Fig. [4f g-i) . In these figures the 
hue shows the polarization eg, red means fully x-polarized, green means fully y- 
polarized and blue is non-polarizing (i.e. equal coupling to both polarizations). 
The intensity of the color shows the Purcell enhancement for that region scaled 
with an enhancement assuming perfect alignment with the field e.g. scaled by 
Fp(x,y) = Fp x (x,y) + Fp y (x,y). From these plots it is clear that for each 
position within the cavity, a preferred dipole polarization exists. Blue areas, 
indicating an equal coupling, only occur where there are nodes in intensity. 
For no point in the cavity is there a region of non-zero intensity which couples 
equally to all polarizations. 

Alternatively, we can calculate the angle of the dipole that gives the maximum 
enhancement. To do this we define the dipole moment as fi = [cos(0); sin(0)], 
where = corresponds to "H" polarization, and 8 = 90° corresponds to "V" 
polarization, and we apply this to Eq. (|14j) . By differentiating this equation 
with respect to 9 and setting it to zero we can find the maximum enhancement 
angle 9 max : 

/j.N 1 , /Im{Gi,2(ro,r ) + G 2 ,i(T , o,r )}\ , . 

ma x(r o ) = - arctan — — — — — — - — (18) 

2 \Jm{Gi,i(r ,ro) + G 2 ,2{ro,r )} J 

where Gij is the matrix element of the GDF on the i-th row and j-th column. 



12 



enhancement 




Figure 4: Spatial map of the Purcell factor for L3 cavities for (a) a horizontal, x, dipole (b) 
a vertical, y, dipole (c) a diagonal, d, dipole (d) an anti-diagonal, a, dipole (e) a right, r, 
dipole and (f) a left, 1, dipole. (g) Degree of polarization for an x and y dipole (h) degree of 
polarization for the d and a dipole and (i) degree of polarization for a r and 1 dipole (note 
that the DOP in (i) is almost zero). 
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Figure [SJa) shows the maximum enhancement angle for the L3-cavity, while 
Fig. [5jb) shows the maximum possible enhancement F max (x,y). The center 
of the cavity shows a maximum enhancement for a dipole at an angle of it/ 2 
which confirms the plots shown in Fig. B|a-f). A few singularities can also 
be seen: these are points where the angle is undefined. One such point is 
indicated by a black arrow. These singularities show points in the PhCC where 
all dipole orientations are enhanced equally, e.g. the eigenvalues of the DGF at 
this point are equal. This would suggest that a spin-photon interface could be 
implemented in L3-cavities. However, if we look at the Purcell enhancement for 
dipoles rotated to the maximum enhancement angles as shown in Fig. [SJb) , we 
observe there is no enhancement at these points, e.g. the eigenvalues are zero. 
QDs positioned at these points are not able to emit into the cavity mode and can 
only emit into leaky modes, this greatly inhibits the spontaneous emission. To 
visualize this we combine Fig.(5Ja-b) to Fig.(S]Jc) where the intensity of the angles 
in Fig. EJa) have been scaled by the enhancement. In this plot the singularities 
have disappeared and no points of orientation independent coupling can be seen. 
This behaviour makes the L3 cavity unsuitable for spin-photon interfaces or the 
entangling photon sources using the biexciton cascade of a QD. 




-k/2 -ti/4 jt/4 7t/2 max max 



Figure 5: (a) Angle of the dipole for which maximum enhancement is achieved, (b) Enhance- 
ment of the dipole when it is oriented at the optimal angle shown in (a).(c) Combination of 
plots (a) and (b), where the angles in (a) have been normalised to intensity using (b). This 
indicates that all areas of degenerate coupling do not show any enhanced emission. 



4.2 Enhancement in HI cavities 

Analysis of the polarization-dependent Purcell enhancement of the dipole mode 
in the HI cavity is somewhat more complex than in the L3 cavity. This is due 
to the presence of two degenerate modes. As we have discussed, the Purcell 
enhancement arises from the QD dipole exciting the mode, which, classically, 
backscatters the light onto the dipole, increasing stimulated emission. Calcu- 
lation of this dipole-mode coupling is relatively easy for an L3 cavity, as only 
one mode is excited. However, for an HI cavity, the dipole couples in general to 
two modes. Thus, whilst the two degenerate modes are orthogonal (as shown 
in Fig. [3]), the dipole may couple the two modes, giving rise to more complex 
backscattering of both modes onto the dipole. 
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For instance, consider the x an d ip modes of the cavity, which were generated 
in Section [3] with a dipole placed in the center, and oriented along x and y, 
respectively (Fig. [3[a-d)). We are able to make this statement because an x 
dipole will only couple to the x component of any mode in the cavity. Examining 
Fig. (3Jc) more closely, we see that the tp mode has a node for the E$(0, 0) = 
component in the center. It will therefore not be excited by an cc-oriented dipole. 
The same goes for a y-oriented dipole: no x mode is excited because, as shown in 
Fig.[3Jb), E*(0, 0) = 0. We can therefore be sure that the two modes generated 
are orthogonal to each other. A dipole oriented along the diagonal, however, 
would excite both the x an d ip modes equally, coupling them together. Thus 
the description of the x an d ^ modes as orthogonal no longer holds when a 
quantum emitter dipole is added: the dipole-mode system should be considered 
as a whole. In fact, as we have discussed, the photonic system is better described 
by a superposition of two orthogonal modes rotated around the "Poincare -like" 
sphere in Fig. Efe). This is what is described by the DGF. 

In general, one should be able to determine which mode will be excited 
by a dipole with an arbitrary orientation at any arbitrary position. As we 
have discussed in Section [3l the HI cavity modes may be expressed as any 
pair of orthogonal modes, arising from a linear superposition of the x an d "0 
modes. These modes can be written in terms of the two basis modes and a 
rotation operator which maps the modes onto two new orthogonal modes on 
the Poincare-like sphere in Fig. 



[ E? } = 



cos £ — sin £ 
sin £ cos £ 



[ E*- E^ ] (19) 



where E*(x, y) and (x, y) are the two new orthogonal basis modes and £ is an 
angle which represents the angle of rotation around the equator of the Poincare- 
like sphere, as shown in Fig. (He). Which basis modes are excited depends on 
both the source position and dipole orientation. As we have seen, when the 
dipole is positioned in the center, rotation of the dipole in the plane by angle 
excites a superposition of the x an d tp modes that is represented by rotation on 
the Poincare-like sphere of £ = 6. This is not the case generally. An arbitrarily 
oriented dipole emitter at an arbitrary position will result in a breaking of the 
symmetry of the system. Thus, rotation of the dipole on the Bloch sphere 
equator does not represent a equal rotation on the Poincare-like sphere, and 
one must find a technique to identify the excited mode. We use Eq. to find 
a rotation angle £ where one of the field components of the modes E^ (r) = 0, 
where p indicates the field component parallel to the orientation of the dipole. 
If -EJ (r) =0, the dipole emitter at r must be at a node of E^ . The orthogonal 
mode of the pair, E*, will give the correct angle of rotation for the system. 

If we define the two eigenmodes in this way, we find that the dipole will show 
no coupling to the E^ mode, coupling only to the E^ mode. The angle £ for 
an arbitrarily polarized source at position r may be expressed as: 

£ p (r) = arctanf^^j (20) 
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The spatial maps £ p (r) for angles of the mode excited by x, y, diagonally and 
anti diagonally oriented sources are shown in Fig. [6] These essentially show 
which type of mode is excited when a dipole of the given orientation (e.g. x) 
is at position tq. One can see here that the mode excited is dependent on 
the QD dipole position. For instance, an x dipole in the center, as shown in 
Fig. [UJa), excites purely the \ mode (£ x = 0), whilst in the "red" region above 
the center the 2-dipole couples purely to the ip mode (£ a = 7r/2). Note that this 
can give some counterintuitive results: in general, one can make no conclusions 
about the orientation of the dipole in the cavity simply by examining the mode 
excited: a vertically polarized dipole emits into the ij) mode if it is in the center, 
but a horizontally polarized dipole also emits into the tp mode if it is off-center. 
Figure E^a-d) also show singularities, i.e. points where the angle is undefined. 
Two such points have been indicated by arrows in Fig. [6ja) . If a dipole is placed 
at these points, no mode can be excited, which would result in inhibition of the 
spontaneous emission. A dipole rotated by £ + n/2 placed at this singularity 
will, however, couple. 




-7i -tc/2 jc/2 n 



Figure 6: Angles of the excited modes for a source placed at a specific location in the 
HI cavity, (a) ^(r) (b) £y(r) (c) §d( r ) (d) §a(f) and (e) Poincare-like sphere indicating the 
angle £ 

Although it is convenient to be able to predict the orientation of the mode 
depending on the dipole position it is not a necessary part of the analysis when 
using DGFs. The DGF for the HI cavity can be calculated using Eq. (fT5|) and is 
independent on the basis chosen for the cavity modes. It is, however, a necessity 
when performing the analysis using Fermi's golden rule. 

The DGF for the HI cavity can be calculated using Eq. (fl5l) using two 
arbitrary orthogonal orientations of the dipole mode. From this, and using 
Eq. (|13p we can plot the Purcell enhancement for any chosen dipole orientation, 
H, as seen in Fig. UJa-f). Figure [Zfg-i) show the DOP, as defined in Eq. (ITT)) , 
for all three bases of dipole orientations in Fig. [7Ja-f ) (rectilinear, diagonal and 
circular). The center of the cavity for all sources show a DOP of zero. This 
suggests that the response of the dipole to the cavity mode is polarization- 
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insensitive. This is in contrast to the L3 cavity, where the DOP is zero only 
for a diagonal/anti-diagonal source in the center of the cavity. Finally, the 
DOP r .i = at all points. This is expected, as no chirality in the cavity structure 
is present. 
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Figure 7: Spatial map of the Purcell factor for HI cavities for (a) Fp x (b) Fp y (c) Fp^ (d) 
Fp a (e) Fp r and (f) Fpi. (g-i) indicate the degrees of polarization for X/Y, D/A and R/L 
dipoles. where the hue shows the degree of polarization and the value (intensity) indicates the 
enhancement 



As with the L3 cavity, we may also calculate the optimal dipole orienta- 
tion for maximum Purcell enhancement using Eq. (|18|) . Figure [SJa) shows the 
optimal angle for the dipole orientation. Using the preferred angles we can now 
continue to plot the maximal enhancement possible in the HI cavity as seen 
in Fig. [5fb), simply by calculating the Purcell factor for the optimal dipole 
orientation in Fig. [5Ja). Here, we note that the Cq v symmetry is preserved. 
Furthermore, the only position with both a degenerate enhancement response 
and a significant Purcell enhancement is in the center. The ring singularity and 
the six singularities between the shifted central holes occur at E-field nodes. 

Next we show the maximum DOP for the preferred angles from Fig. EJa) in 
Fig. [HJc). In this plot, the singularities that we mentioned in Fig. Eta), show 
up again in the form of a DOP of zero, meaning there is no preferred angle of 
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0.25 O.S 0.75 



Figure 8: (a) Angle of the dipole for which maximum enhancement is achieved, (b) Enhance- 
ment of the dipole when it is oriented at the angle shown in (a), (c) maximal degree of 
polarization and (d) a combination of (a) and (b) where the intensity of the plot is scaled by 
the enhancement. 



emission. From this plot it is easy to determine where a QD dipole would couple 
strongly to one polarization only. 

Finally we combine Fig. EJa-b) in one plot (as in Fig. EJc)) and can be seen 
in FigJHJd). From this plot it becomes clear that the dipole's response to the 
cavity is polarization insensitive if it is positioned perfectly in the center of the 
cavity. 

When the QD is not perfectly positioned in the center of the cavity the spin- 
photon interface will not have unit fidelity. The actual fidelity of the interface 
is highly dependent on the physical underlying principle and the read in and 
read out channels. In order to give some sort of length scale, we give the ratio 
of the minimum (FP min ) and maximum (FP max ) enhancement, by calculating 
the Purccll enhancement at the optimum angle and perpendicular to it, from 
the values in Fig.[8ta). This is an indication of the minimum fidelity that should 
be obtainable. We define this minimum fidelity as: 




□ 0.25 0.5 0.75 



Figure 9: (a) Contour plot of the minimum fidelity of the spin-photon interface, (b) Maximum 
fidelity when the QD is moved along the x-axis from the center of the cavity. 

Figure Eta) shows a contour map of the fidelity measure from Eq. (|2Tj) . A 
few spots can be seen where the fidelity is close to unity, however, as mentioned 
earlier, only the central spot shows an enhancement and is the only suitable 
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QD location. In order to give a required precision for the QD location, we 
plot the fidelity against the QD position along a line in the x-direction as can 
be seen in Fig. HJb). If a fidelity of 95% is required, the QD needs to be 
positioned within 0.09a of the center. Assuming a typical resonant wavelength 
of 900nm, a lattice constant of 260nm is required for our HI design. This 
means that the QD needs to be positioned within approximately 20nm of the 
center of the cavity. Note that this positioning accuracy is more stringent than 
the ~ 60nm estimated to be needed for strong coupling [35]. This criterion 
pushes the limits of what has presently been demonstrated for site controlled 
QDs [2], but methods using QD identification and PhC positioning techniques 
have shown that the combined QD-PhC lithography positioning accuracies of 
~ 30nm approach these requirements [35] . 

5 Conclusion 

In conclusion, we have demonstrated a method to determine the suitability of 
photonic structures for a QD spin-photon interface. We addressed the problem 
from the specific point of view of QD-PhC structures for quantum information, 
by considering firstly in which cavity designs one is able to to create inter- 
mediate quantum superposition states, and showed that a doubly-degenerate 
cavity mode is required. We then consider the influence of the QD position 
on the ability to excite a superposition state. We demonstrate that even in 
the degenerate HI cavity, QD position is as crucial to the superposition state 
formed as the dipole orientation. In particular, to allow accurate transmission 
of a QD spin superposition state to a cavity superposition state, the QD must 
be accurately placed to within ~ 20nm. 

We note that the methods described in this paper may be applied to a 
quantum emitter placed in any type of lossless photonic cavity environment 
that modifies the LDOS, in order to calculate polarization-dependent Purcell 
factors for any emitter position and orientation. In particular, more complex 
cavity designs (such as coupled cavity, or cavity waveguide designs) may be 
analyzed using simple FDTD simulations and these techniques for suitability 
for generating and transmitting superposition states. 
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